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Abstract 

Porteresia coarctata is a wild relative of rice with capabilityof high salinity and submergence tolerance. The 
transcriptome analyses of Porteresia can lead to the identification of candidate genes involved in salinity and 
submergence tolerance. We sequenced the transcriptome of Porteresia under different conditions using 
lllumina platform and generated about 375 million high-quality reads. After optimized assembly, a total of 
1 52 367 unique transcript sequences with average length of 794 bp were obtained. Many of these sequences 
might represent fragmented transcripts. Functional annotation revealed the presence of genes involved in 
diverse cellular processes and 2749 transcription factor (TF)-encoding genes in Porteresia. The differential 
gene expression analyses identified a total of 15 158 genes involved in salinity and/or submergence 
response(s). The stress-responsive members of different TF families, including MYB, bHLH, AP2-EREBP, 
WRKY, bZIPand NAC, were identified. We also revealed key metabolic pathways, including amino acid biosyn- 
thesis, hormone biosynthesis, secondary metabolite biosynthesis, carbohydrate metabolism and cell wall 
structures, involved in stress tolerance in Porteresia. The transcriptome analyses of Porteresia are expected 
to highlight genes/pathways involved in salinity and submergence tolerance of this halophyte species. The 
data can serve as a resource for unravelling the underlying mechanism and devising strategies to engineer sal- 
inity and submergence tolerance in rice. 
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1 . Introduction 

High salinity is a severe and increasing constraint on 
the global crop production. Most of our major crops 
are highly susceptible to salinity, and there is an urgent 
need to address this problem. Salinity tolerance is a 
complex multigenic trait, which depends on a range of 
physiological processes. The salinity stress response has 
been studied in model crop plant rice,' which repre- 
sents a glycophyte plant and does not possess the cap- 
ability to tolerate high salinity. The existence of wild 
halophytic species, which can grow not only in saline 



but also in submergence prone habitats, has been 
reported. ^'^ These wild species may prove to be an im- 
portant gene reservoir for improving salinity tolerance 
in crop plants via genetic engineering approaches.^"^ 

Porteresia coarctata {Oryza coarctata) is a wild halo- 
phyte, which shows considerable adaptation to high 
salinity (20-40 dSm~') and can also tolerate com- 
plete submergence in saline water for a long period.^ 
Porteresia grows in abundance near coastal areas, 
where the vegetation is inundated twice a day with 
highly saline water. Although many salinity- and sub- 
mergence-tolerant rice cultivars have also been 
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identified,' ^theyseemto be much lesstolerantthan 
Porteresia. It is a close wild relative of rice and can serve as 
a n i m porta nt sou rce of genes related to sa I i n ity a nd su b- 
mergence tolerance.^'' ^ Some novel insights into the 
molecular mechanisms underlying salinity tolerance 
have been obtained from transcriptome and genome 
analyses of close halophyte relatives of Arabidopsis, 
Thellungiella halophila and T. parvula}'^'^^ A few 
studies have demonstrated that engineering of some 
genes from Porteresia, such as those involved in inositol 
metabolic pathway, can confer salinity tolerance in cul- 
tivated rice.' ^ However, only a few genes have been 
isolated from Porteresia so far. Furthermore, the effect 
of submergence stress and its combined effect with 
salinity have hardly been analysed.^ The discovery of 
more genes and metabolic pathways in wild species 
like Porteresia can provide a better strategy to impart 
salinity tolerance in cultivated rice. 

The availability of next-generation sequencing 
technologies provides a high-throughput, rapid and cost- 
effective means to sequence and characterize the tran- 
scriptome of non-model and wild species as well.'^~^' 
In this study, we sequenced the transcriptome of 
Porteresia underdifferent conditions using llluminatech- 
nology. The transcriptome assembly was generated 
to reveal the first picture of its gene content. We also ana- 
lysed gene expression under salinity and submergence 
conditions and identified the genes and metabolic 
pathways involved in these stress responses. The data 
presented here will provide a very useful resource to 
engineer salinity and submergence tolerance stress in 
crop plant rice. 

2. Materials and methods 

2.7. Plant material 

Porteresia coarctata (Roxb) Tateoka plants were 
collected from saline river banks (Matia) of the 
Sunderbans in the coastal region of Bay of Bengal, 
India. The plants were brought to the laboratory along 
with soil and water. The plants were washed carefully 
to remove soil from the leaves and roots. The plants 
were given different treatments, including control (RO 
water, control), low salinity (450 mM NaCi, salt450), 
high salinity (700 mM NaCi solution, salt700), submer- 
gence (fully submerged in RO water, submergence) and 
submergence in low salinity (fully submerged in 
450 mM NaCi, salt + submergence). The plants were 
harvested in three biological replicates 1 2 h after the 
treatment. 

2.2. RNA isolation and sequencing 

Total RNA was isolated from each tissue sample 
(control,salt450,salt700,submergence and salt + sub- 
mergence) using TRI reagent (Sigma Life Science, USA). 



The quantity and quality of RNA samples were checked 
using Nanodrop (Thermo Fisher Scientific) and Agilent 
Bioanalyzer (Agilent Technologies, Singapore). Equal 
amounts of total RNA from three biological replicates 
were pooled for all the five treated samples before 
library preparation. Sequencing of all the five samples 
was performed by the commercial service provider, 
Genotypic Technology Pvt. Ltd, Bangalore, India, each 
in one lane using an llluminaGenome Analyzer II to gen- 
erate 1 00 base long paired-end reads. The raw data 
obtained were filtered through the standard lllumina 
pipeline. The filtered Fastq files were further subjected 
to more stringent quality control using the NGS QC 
Toolkit (v2.3) to remove the low-quality reads and 
reads containing adaptor/primer contamination.^^ 
The high-quality filtered reads were used for further 
downstream processing. The entire sequence data gen- 
erated in this study are available in the Gene Expression 
Omnibus database under the accession number 
GSE4491 3. The non-redundant lllumina reads were fil- 
tered using the CDHIT program (http://weizhong-lab. 
ucsd.edu/cd-hit/). 

2.3. Transcriptome assembly 

Denoi^o transcriptome assembly was performed using 
various tools, including Velvet (vl.2.01). Oases 
(vO.2.04), ABySS (vl .2.6), Trinity (vr201 2-05-1 8) and 
CLC Genomics Workbench (v4.7.2). Various parameters, 
including k-mer length, insert length and expected 
coverage, were optimized to obtain best assembly as 
described previously.^^ Reference- based assembly 
was performed using the CLC Genomics Workbench 
using default parameters, and assembled consensus 
sequences were exported into a fasta file. The annotated 
cDNA sequences of japonica rice (Os-Nipponbare- 
Reference-IRGSP-1 .0)^'' available at Rice Genome 
Annotation Project (MSUv7.0; http://rice.plantbiology. 
msu.edu/) were used as a reference. Merged assembly 
was performed using the TGICL program (v2.0), with 
minimum overlap length of 40 and minimum percent- 
age identity of 90 for overlaps. The quality of assemblies 
generated was assessed using the perl scripts given in the 
NGS QC Toolkit. 

2.4. Functional annotation 

BLASTX searches against rice proteome, UniRef90, 
UniRefl 00 and NCBI non-redundant nucleotide data- 
bases, were performed to assign a putative function to 
each Porteresia transcript. GOSlim terms to each tran- 
script were assigned on the basisof their bestsignificant 
match (E-value < 1 e- 1 0) rice protein to provide a 
broad overview of their function. An E-value cut-off of 
le-10 was used to identify orthologs between 
Porteresia and rice. The transcription factor (TF) encod- 
ing Porteresia transcripts were identified based on the 
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hidden Markov model profile search of conserved 
domain presentineachTFfamilyasdescribed previous- 
ly.^^ The statistical significance of the difference in the 
number of members between Porteresia and rice in 
each TF family was calculated via Fisher's exact test. 

2.5. Read mapping and gene expression analysis 

To estimate the gene expression, all the high-quality 
reads from each condition were aligned on the tran- 
scriptome assembly using the RNA-Seq Analysis utility 
of CLC Genomics Workbench. A maximum of two mis- 
matches were allowed for the alignments. The read 
countsforeach tissue sample were normalized by calcu- 
lating the read per kilobase per million for each tran- 
script. The DESeq software (vl.l 0.1)^^ was used for 
differential gene expression analysis. DESeq provides a 
simple and computationally faster method for quanti- 
tative measurement of gene expression based on the 
negative binomial distribution with variance and 
mean linked by local regression. We calculated the 
size factor for each sample for normalization of read 
count data using DESeq. A P-value cut-off of <0.05 
and at least two-fold change was used to identify differ- 
entially expressed genes. Heatmaps showing expression 
profiles (log2 fold change) were generated using the 
MultiExperiment Viewer (MeV, v4.8). /<-means cluster- 
ing was performed using Euclidean distance matrix 
with 1 000 iterations using MeV. Hierarchical clustering 
was performed using Euclidean distance matrix with a 
complete linkage rule using MeV. 

2.6. Gene ontology and pathway enrichment analysis 
The best Arabidopsis hit corresponding to each 

Porteresia transcript was identified using BLAST search 
to study gene ontology (GO) enrichment. GO enrich- 
ment of different sets of differentially expressed genes 
was performed using the BiNGO tool.^^ The metabolic 
pathway data available in the RiceCyc database of 
Gramene^^ were analysed to identify the enriched 
metabolic pathways in various gene sets. Both the GO 
and pathway enrichment analyses were performed at 
P-value cut-off of <0.05 after applying Benjamini- 
Hochberg correction. 

2.7. Real-time PGR analysis 

For real-time PGR analysis, gene-specific primers 
were designed using the Primer Express (v3.0) software 
(Applied Biosystems, Santa Clara, CA, USA). The primer 
sequences used in this study are listed in 
Supplementary Table SI. Real-time PCRs were per- 
formed as described. The expression of most suitable 
internal control gene for rice, UBQS,^^ in each sample 
was used for normalization. At least three independent 
biological replicates with three technical replicates of 



each biological replicates for each tissue sample were 
used for the analysis. The correlation between expres- 
sion profiles of selected genes obtained from real- 
time PGR and RNA-seq data analyses was determined 
using MS Excel. 

3. Results 

3.1 . Transcriptome sequencing 

To generate the transcriptome of Porteresia, we col- 
lected tissue samples from plants kept in water 
(control) and those subjected to salinity and sub- 
merged conditions. For salinity stress, we kept the 
plants in 450 mM (salt450) and 700 mM (salt700) 
sodium chloride solutions and for submergence stress, 
plants were kept completely submerged in water (sub- 
mergence) and 450 mM sodium chloride solution 
(salt + submergence) and collected the tissues after 
1 2 h of treatment. Porteresia is a halophyte, which 
grows naturally in coastal region, where plants experi- 
ence lunar tide and are exposed to submergence in 
saline sea water every 1 2 h. Therefore, to obtain 
maximum representation of genes involved in signal- 
ling and adaptation response under salinity and sub- 
mergence stress, we selected time point of 1 2 h for 
the analysis in our study. The RNA-seq libraries were 
constructed and sequenced using an lllumina sequen- 
cing platform. We generated a total of more than 422 
million reads from above five samples. After quality fil- 
tering, about 375 million high-quality reads (ranging 
from about 65 to 90 million reads for each sample) 
were obtained (Table 1). The average Phred quality 
score was at least 30 at each base position for filtered 
reads, indicating high quality of filtered reads 
(Supplementary Fig. SI ). 

3.2. Transcriptome assembly optimization 

To generate the optimal assembly of the Porteresia 
transcriptome, first, we performed de novo assembly 
using all the high-quality reads (374 759 624) and 
non-redundant reads (1 94 31 6 505) after removing 
duplicate reads from all the five samples. We used 
Velvet, Oases, ABySS, Trinity and CLC Genomics 
Workbench softwares to perform de novo assembly. 
The de novo assembly was generated at different 
k-mer lengths using these softwares. From each soft- 
ware, the best assembly output was selected based on 
the statistics (average and N50 contig lengths) and 
similarity with rice proteome (Supplementary Table 
S2). Among the various softwares. Oases generated 
the best assembly with largest average and N50 
lengths at /<:-mer length of 59 for both datasets 
(Supplementary Table S2). However, only 34% and 
29% rice proteins showed a significant hit with the 
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Tot3l reads 


High-quality reads 


Mapped reads^ 


Mapped reads (%) 


Control 


73 200 1 60 


64 799 692 


56 638 1 1 3 


87.40 


Salt450 


78 71 7 1 62 


68 275 254 


59 81 8 916 


87.61 


SaltZOO 


73 664 082 


66 857 1 02 


58431 827 


87.40 


Submerged 


1 01 51 7 546 


89 800 570 


78 997 383 


87.97 


Salt450 + submerged 


95 241 740 


85 027 006 


74 867 541 


88.05 


Total 


422 340 690 


374 759 624 


328 753 780 


87.72 



Number of high-quality reads mapped using CLC Genomics workbench. 

Table 2. Assembly statistics of Porteresia transcriptome 



assemblies generated from whole and non-redundant 
datasets, respectively. The assembly output from other 
softwares also exhibited very low similarity to rice pro- 
teins. As Porteresia is a wild relative of rice, higher simi- 
larity to rice proteins was expected. These results 
suggested that the output of de novo assembly was not 
optimal. 

Therefore, we followed a combined approach of both 
reference- based and de novo assemblies to generate the 
transcriptome of Porteresia. In the first step, we per- 
formed reference-based assembly using the CLC 
Genomics Workbench. We used all the reads from 
each sample to generate a consensus assembly 
(Supplementary Table S3). The average and N50 
lengths of reference assembly were much higher 
when compared with de novo assemblies. The consen- 
sus assembly from all the five samples generated a 
total of 36 966 contigs with the average length of 
1 265 bp and N50 length of 1 761 bp (Supplementary 
Table S3). This assembly utilized ~81% of the total 
reads. To identify the novel genes in Porteresia, as a 
second step, we performed de novo assembly of 
un used reads i n the first step of reference-based assem- 
bly. We used different softwares (CLC, Trinity, ABySS, 
Velvet and Oases) and parameters to select the best as- 
sembly. Oases generated the best assembly with largest 
average (544 bp) and N50 (719 bp) contig lengths. 
However, the number of consensus sequences 
(1 61 1 60) generated was quite large. As a last step, 
we generated a merged assembly of the contigs gener- 
ated in primary assemblies of two steps via TGICL 
program. Our earlier studies have shown that TGICL 
generates optimal merged assembly from the contigs 
of primary assemblies.^^'^° This assembly resulted in a 
totalof 1 52 367 (39 203 contigsand 113 1 64 single- 
tons) unique sequences (referred as contigs) of at least 
200 bp in length (Table 2). The average and N50 
lengths of contigs were 794 and 1 249 bp, respectively. 
The mapping of all the high-quality reads suggested 
that > 9 5% of reads were used in the final merged as- 
sembly and average read depth of each contig was 
2352. About 55% of rice proteins were represented in 
this assembly and >23% of rice proteins exhibited 
> 70% coverage (Table 2). 



Total reads 


374 759 624 


Total transcripts 


1 52 367 


Size (Mb) 


1 20.99 


Minimum transcript length (bp) 


200 


Maximum transcript length (bp) 


14945 


Average transcript length (bp) 


794.07 


N50 length (bp) 


1 249 


Total reads used for assembly 


95.62% 


Average read depth 


2352 


Transcripts with significant hits" 


48 048 (31 .5%) 


Rice proteins with significant hit'' 


21 499 (55.06%) 


Rice proteins with > 70% coverage"^ 


91 14 (23.34%) 


^Number of transcripts showing a si^ 


;nificant hit (£-value < 


1 e- 5) with rice proteins. 




''Number of rice proteins showing a si; 


gnificant hit (£-value < 


1 e - 5) with Porteresia transcripts. 




"^Number of rice proteins that showed > 70% coverage. 



3.3. Functional annotation o/Porteresia transcriptome 
The unique set of transcripts generated above were 
designated as P. coarctata tentative consensus (PcTC) 
transcripts and assigned unique identifier number 
from PcTCOOOOOl to PcTCl 52367. The whole 
transcriptome sequence is available at Porteresia 
Transcriptome Sequence web page (http://nipgr.res. 
in/mjain.html?Page=porteresia). The total transcrip- 
tome assembly size was ~1 21 Mb. About 48% of all 
the transcripts were longer than 500 bp and >24% 
transcripts were larger than 1000 bp (Supplementary 
Fig. S2). Most of the smaller transcripts were repre- 
sented by the singletons (Supplementary Fig. S2). 
Many of these are expected to merge into larger tran- 
scripts, once more sequence data of larger read length 
become available. The average GC content and distribu- 
tion of Porteres/fl transcripts were significantly different 
than that of rice (Supplementary Fig. S3), indicating 
substantial differences between rice and Porteresia 
genomes. Furthermore, to explore the putative func- 
tion, we assigned gene description to each Porteresia 
transcript corresponding to the ortholog rice proteins 
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and proteins in other databases. A gene description term 
could be assigned to ~ 3 6% of Porteres/fl transcripts. Tfie 
lower percentage of annotated transcripts was because 
most of the smaller transcripts generated in the de novo 
assembly did not show significant similarity with 
known proteins in various databases and presumably 
represent species-specific genes. 

Based on their similarity with rice genes, Porteresia 
transcripts were assigned GOSlim terms under biologic- 
al process, molecular function and cellularcomponent 
categories. Amongthe biological process terms, biosyn- 
thetic process (20%) was most represented followed by 
nucleicacid metabolic process (1 5.3%) and response to 
stress (14.5%) (Fig. 1A). A high representation of re- 
sponse to stress category might contribute to the high 
salinity and submergence tolerance phenotype of 
Porteresia. The largest number of transcripts belonged 
to catalytic activity (1 6.1 %), protein binding (1 5.1 %) 
and hydrolase activity (1 2.3%) among the molecular 
function GOSlim terms. Amongthecellularcomponent 
terms, the transcripts related to membrane (18.5%) 
and plastid (1 7%) were the largest in number 
(Fig. 1 A). We further analysed the Porteresia transcrip- 
tome to identify the TF-encoding genes. A total of 
2 749 transcripts encoding for TFs were identified. 
Amongthe 85 families represented in these transcripts, 
theMYB-domain (1 97)family TFs were mostabundant 
followed by bHLH (1 50)- and NAG (1 47)-domain TFs 
(Fig. 1 B). Although the total numberof TFs was slightly 
higher in Porteresia (2 749) than in rice (2452), we 
found significant differences in their distribution in 
certain families. Notably, the number of TFs predicted 
in Porteresia for bromodomain, zinc-finger CCHC and 
SNF2 families were significantly higher than rice. 
Likewise, the number of predicted TFs in AP2-EREBP, 
GRAS, MADS, TCP, TRAF and ZF-HD families were signifi- 
cantly lesser in Porteresia than in rice (Fig. 1 B). These 
results suggest significant differences in the transcrip- 
tional machinery of the twospeciesand might contrib- 
ute to stress adaptation phenotype in Porteresia. 

3.4. Differential gene expression analysis 

In addition to the transcriptome construction, RNA- 
seq provides an opportunity to measure gene expres- 
sion at the whole-genome level.^''^^ To investigate 
the differential gene expression under different condi- 
tions, we mapped the reads from each sample on 
the Porteresia transcriptome. A total of 87-88% of the 
high-quality lllumina reads were mapped to the 
Porteresia transcripts (referred as genes hereafter) 
from each sample (Table 1 ). Differential gene expres- 
sion analysis was performed using the DESeq software 
to identify the genes with altered expression under dif- 
ferent conditions. Atotal of 1 5 1 58 genes exhibited dif- 
ferential expression under at least one stress condition 
analysed (Supplementary Table S4). The largest 



number of genes was differentially expressed in the 
plants subjected to submergence stress (5864) fol- 
lowed by both salinity and submergence stresses to- 
gether (5787) (Fig. 2A). A total of 4841 and 4243 
genes exhibited differential expression under low salin- 
ity (salt450) and high salinity (salt700) stress condi- 
tions, respectively. Under salinity stress, a larger 
number of genes were up-regulated (3020 under 
salt450 and 2553 under salt700) when compared 
with down-regulated genes (1 821 under salt450 and 
1 690 under salt700). However, in contrast, the larger 
number of genes was down-regulated under submer- 
gence stress alone (3712) and salinity and submer- 
gence stresses together (3952) when compared with 
up-regulated genes (2152 for submergence and 
1 835 for salt + submergence) (Fig. 2A). 

Furthermore, we identified genes showingoverlapping 
or specific response to different stress conditions as 
depicted in Fig. 2B. Interestingly, a larger fraction of 
genes (69.5%) exhibited response to a specific stress con- 
dition. Other 4625 (30.5%) genes were responsive to 
more than one stress conditions. The largest number of 
genes (2 555) wascommonly regulated by submergence 
and salt + submergence stresses. Notably, most of them 
(99.6%) exhibited similar response (681 up-regulated 
and 1865 down-regulated) under both conditions. 
However, among the 790 genes commonly regulated 
by both high salinity and submergence stresses, most 
(7 8.2%) exhibited opposite responses (473 up-regulated 
under high salinity but down-regulated under submer- 
gence stress and 145 down-regulated under high 
salinity but up-regulated under submergence stress). 
Interestingly, however, 84%of the 757 genes common- 
ly regulated by low salinity and submergence stresses 
exhibited similar response. Only 957 genes exhibited 
differential expression under both low and high salinity 
conditions, but most of them (94.6%) were regulated in 
thesamedirection.Averysmallfractionof genes(l 1 6) 
were differentially regulated under all the stress condi- 
tions analysed. Furthermore, we performed /<-means 
clustering analysis and generated 1 0 clusters with 
different gene expression patterns across different 
stress conditions analysed (Supplementary Fig. S4). 
Among the clusters generated, we identified four 
groups of genes with distinct gene expression patterns 
(Fig. 3). Group I genes were up-regulated under salinity 
and/orsalt + submergence stress. Group II genes were 
down-regulated under submergence and salt + sub- 
mergence stress. Group III genes were up-regulated 
under submergence and salt + submergence stress. 
Group IV genes were up-regulated under salinity stress. 

3.5. Stress response of transcription factor genes 

TFs are the regulatory proteins, which play a key role 
in stress responses. Remarkably, we detected the 
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Figure 1 . Functional annotation of Porteresia transcripts. (A) GOSlim term assignment to the Porteresia transcripts in different categories of 
biological process, molecular function and cellular component. (B) The number of transcripts representing different TF families in 
Porteresia and rice transcriptomes are shown. Fisher's exact test was used to evaluate the significant difference in the number of 
members representing TF families in the two species. The TF families with statistically significant difference in the number of members 
are marked with single (P < 0.05) or double (P < 0.01 ) asterisks. 



differential expression of at least 329 (~12%) TF- 
encoding genes under different stress conditions in 
Porteresia. A heatmap showing the differential 



expression of genes belonging to different TF families 
isshown in Fig. 4. At least one memberofSSTF families 
exhibited differential expression. Among the various 
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Salt450 Salt700 Submergence Salt + 

submergence 




Figure 2. Differential gene expression under different stress 
conditions. (A) The number of up- and down-regulated genes 
under different stress conditions (salt450, saltZOO, submergence 
and salt + submergence) is shown in the bar graph. The total 
number of genes differentially expressed under each condition is 
given on the top of each ban (B) Venn diagram showing the 
number of genes differentially expressed under specific or 
different stress conditions. 

families, the largest number (32) of NAC TF family 
genes was differentially expressed. The members of 
other TF families, including MYB/MYB-related (2 5), 
bZIP (25), bHLH (23), HSF (1 8) and WRKY (1 7), were 
also significantly represented in the differentially 
expressed genes. A total of 96, 1 18, 1 32 and 1 27 TF 
genes exhibited differential expression under low salin- 
ity, high salinity,submergence and salt + submergence 
stress conditions, respectively. Among the 1 1 8 TF genes 
differentially expressed under high salinity stress, 88 
were up-regulated and 30 were down-regulated. The 
members of NAC, MYB/MYB-related and WRKYfamilies 
were the most abundant among these genes. In con- 
trast, among the 1 32 TF genes differentially expressed 
under submergence stress, 44 were up-regulated and 



88 were down-regulated. The members of bHLH, bZIP, 
HSF, NAC and MYB/MYB-related families were most 
abundant among these genes. In addition, several 
zinc-finger TFs of different classes were also regulated 
under different stress conditions. The members of 
these TF families have been found to be regulated 
under salinity and submergence stresses in Arabidopsis, 
rice and halophytes.^^ Altogether, these data 
suggest the extensive regulation of transcription under 
salinity and submergence stress. 

3.6. Gene ontology and metabolic pathway enrichment 
analysis of stress-responsive genes 
To investigate the major functional categories repre- 
sented in the stress-responsive genes identified above, 
we performed GO enrichment analysis. Most evidently, 
besides response to stress, we found the GO terms asso- 
ciated with various metabolic processes, such assecond- 
ary metabolite (oxoacid, ketone, amino acid, oxylipin, 
flavonoid, phenylpropanoid and chorismate) biosyn- 
thesisand hormone (jasmonicacid,ethyleneandabsci- 
sic acid, ABA) biosynthesis, were most enriched in the 
genes up-regulated under different stress conditions 
analysed (Fig. 5, Supplementary Fig. S5). In addition, 
the GO terms, photosynthesis, photorespiration, rRNA 
processing and transport were also found enriched in 
up-regulated genes. However, in the down-regulated 
genes under salinity stress, the GO terms associated 
with polysaccharide (glucan, carbohydrate and starch) 
metabolism were most significantly represented 
(Supplementary FigsS5 and S6).The genes with function 
in protein turn over (proteolysis, peptidase activity and 
proteasome assembly) were highly enriched under low 
salinity. Such genes have been found responsible for sal- 
inity stress adaptation in halophyte,T /ifl/op/7;7fl.^^ Under 
submergence alone or along with salinity stress, the GO 
terms response to various stresses, response to ABA, gen- 
eration of precursor metabolites and energy as well as 
rRNA metabolism were most enriched. The finding of 
up-regulation of genes associated with energy conserva- 
tion via metabolism rearrangement to enhance ATP pro- 
duction and/or nullifying deleterious effects of oxygen 
radicals, reactive oxygen species and other toxic com- 
pounds under submergence has been reported earlier 
also.^^'^^ In addition, the GO terms, such as response to 
gibberellin stimulus, response to hydrogen per oxide 
and transport, were represented in the down-regulated 
genes under submergence stress (Supplementary Figs 
S5 and S6). The GO term, regulation of transcription, 
was also found enriched in the stress-responsive genes, 
which is consistent with differential expression of the 
large number of TF genes as described in previous 
section. Notably, we noticed the significant enrichment 
of cellularcomponent GO term, cell wall, in the differen- 
tially expressed genes (Supplementary Fig. S7). The 
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Figure 3. /<-means clustering of expression profiles of genes differentially expressed under different stress conditions. The clustering was 
performed on log2 fold change for each gene under different stress conditions when compared with under control condition. The genes 
exhibiting a similar pattern of expression under different stress conditions [salt450, saltZOO, submergence (sub) and salt + 
submergence (salt + sub)] have been grouped together in four distinct groups. The scale at the bottom shows log2 fold change. A colour 
version is available online. 



changes in transcript levels of a number of genes encod- 
ing cell wall-related proteins have been observed in 
halophyte, T. halophila, in response to abiotic stress con- 
ditions.^^'^^ Overall, the different groups of genes with 
distinct expression patterns were also found to be asso- 
ciated with similar GO terms (Supplementary Fig. S8). 
For instance, the GO terms, response to stress and 
hormones were enriched in Group II genes, whereas 



photosynthesis and secondary metabolite biosynthesis 
terms were most significant in Group III genes. Group IV 
genes were found enriched in tryptophan biosynthesis, 
cellular amino acid derivatives biosynthesis, response to 
jasmonic acid and cell death-related GO terms. 

Many of the stress-responsive genes identified above 
encode for the enzymes involved in several metabolic 
pathways, which regulate diverse cellular processes. 
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Figure 4. Differential expression of TF-encoding genes under 
different conditions. Various TF families showing differential 
expression under different stress conditions [salt450, saltZOO, 
submergence (sub) and salt + submergence (salt + sub)] are 
given on the right side. The scale at the bottom represents log2 
fold change. A colour version is available online. 



Therefore, we investigated the metabolic pathways 
significantly enriched in stress-responsive genes by 
integrating the stress-responsive gene set with meta- 
bolic pathway data available in the RiceCyc database. 
The pathways involved in amino acids biosynthesis, 
secondary metabolite biosynthesis, hormone biosyn- 
thesis, carbohydrate metabolism and cell wall structures 
were significantly overre presented among stress- 
responsive genes (Fig. 6), which are in well agreement 
with GO enrichment analysis. A list of pathways affected 
under stress conditions are given in Supplementary 
Table S5. Most evidently, the genes involved in the bio- 
synthesis of various amino acids, including tryptophan, 
valine, isoleucine, glycine, glutamate and serine, were 
significantly represented in stress-responsive genes. 
The pathways for various stress-related hormone 
biosynthesis, such as jasmonic acid, ethylene and ABA 



biosynthesis, were also enriched. The genes involved in 
jasmonic acid biosynthesis were up-regulated under 
salinity stress, whereas genes for ABA biosynthesis 
were up-regulated under salinity stress, but down- 
regulated under submergence stress (Fig. 6). The path- 
ways involved in carbohydrate metabolism, such as 
starch degradation, trehalose biosynthesis and anaer- 
obic respiration, were also overrepresented among the 
stress-responsive genes. The genes involved in anaerobic 
respiration were significantly up-regulated under sub- 
mergence stress alone and along with salinity stress. The 
genes involved in anaerobic metabolism have been 
found regulated in response to submergence, flooding 
and O2 deprivation across the plant species too.^^' 
Notably, the genes involved in cell wall structure bio- 
synthesis, such as cellulose and suberin biosynthesis, 
were also significantly enriched. The pathways involved 
in serotonin amides biosynthesis, phenylpropanoid 
biosynthesis, phospholipase activity and fatty acid 
derivative biosynthesis were also found enriched. A 
similarset of metabolic pathways weresignificantly repre- 
sented in the different groups of genes with distinct ex- 
pression patterns (Supplementary Table S6). 

Interestingly, we detected the differential expression 
of several genes in Porteresia implicated in SUB1A- 
mediated submergence tolerance response."^" In 
addition to the genes involved in ethylene response, 
carbohydrate metabolism, photorespiration and 
anaerobic respiration, we found putative gibberellin 
biosynthesis genes regulated under submergence 
stress, wh ich med iates shoot elongation. The expression 
of Alcohol Dehydrogenase, a marker gene for anoxia 
stress, was also up-regulated under submergence 
stress, consistent with previous findings.^^'"^"""^^ We 
detected differential expression of several putative 
expansin genes as well under submergence and/or 
salt + submergence stress. The expansins have been 
shown to play a major role in cell expansion under 

anoxic condition, prevalent under submergence 
stress.^'^'^^-'^^ 



3 . 7. Real-time PCR validation of differential gene 
expression 

We validated the results of differential gene expres- 
sion analysis obtained from RNA-seq data by real-time 
RT-PCRanalysis.Theexpressionofatleast 1 0 randomly 
selected genes was analysed in the three biological 
replicates of five tissue samples representing different 
conditions. The analysis revealed similar expression 
pattern of all the selected genes in real-time PCR ana- 
lysis as observed from RNA-seq data. The statistical 
analysis also showed very good correspondence 
(correlation coefficient of 0.792) among the results 
of real-time PCR and RNA-seq data analyses as shown 
in Fig. 7. 



78 



Transcriptome Analyses in Wild Halophyte Porteresia 



[Vol. 21, 





Figure 5. GO enrichment under different stress conditions. Significantly, enriched GO categories in the up-regulated genes under high salinity 
(A) and submergence (B) stress are shown. The genes were analysed using BiNGO and the biological process terms showing significant 
enrichment are shown. Node size is proportional to the number of genes in each category and shades represent the significance level 
(white— no significant difference; scale, P = 0.05 to P < 0.0000005). A colour version is available online. 
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Figure 6. Regulation of metabolic pathways during salinity and submergence stress conditions. The metabolic pathways enriched in 
differentially expressed genes are shown. Heatmaps showing the expression profiles of genes involved in these pathways under different 
stress conditions [salt450, saltZOO, submergence (sub) and salt + submergence (salt + sub)] are also shown. The scale at the bottom 
represents log2 fold change. A colour version is available online. 



4. Discussion 

Over the past many years, considerable efforts have 
been made to understand the molecular mechanism 
of salinity tolerance/adaptation in plants. Several 
studies have reported the molecular responses of 



plants to salinity stress.^ However, only a few 
studies have been performed to analyse the submer- 
gence response.^^'^^'"^" In the recent years, study of 
wild halophyte relatives has gained attention as models 
to investigate the mechanism of salinity and submer- 
gence tolerance due to their extreme tolerance/ 
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Figure 7. Validation of differential gene expression results obtained 
by RNA-seq. The correlation of gene expression results obtained 
from real-time PGR analysis and RNA-seq for 1 0 selected genes is 
shown. 

adaptability to these conditions.^'^ However, most of 
such studies are confined to close relatives of 
Arabidopsis, T. halophila and T. parvula}'^'^^ Recently, 
the transcriptome of a salinity-tolerant wild tomato 
species has also been reported.'*'^ However, considering 
the complexity of salinity response and fundamental dif- 
ferences in the physiology of monocots and dicots, these 
dicot halophytesare not ideal modelsto understand the 
mechanism of salinity tolerance in monocot crop plants. 
Recently, the importance o\P. coarctata,\N\\6 halophyte 
relative rice, has been recognized as a model to study 
stress tolerance in monocots and important source of 
genes for the improvement of rice via biotechnological 
approaches.^ 

In the present study, we sequenced the transcriptome 
of Porteresia under different stress conditions using 
lllumina technology to explore its gene complement. 
Unexpectedly, Porteres/a transcripts exhibited low simi- 
larity with rice proteins, suggesting significant differ- 
ences between Porteresia and rice genomes. The 
similarity of Porteresia genes with rice was lower than 
that of other sequenced wild rice species.'*^''*^ A sub- 
stantial number of genes appeared specific to 
Porteresia because they do not have an ortholog in 
rice and other organisms as predicted based on the se- 
quence similarity. Furthermore, we observed expansion 
or contraction of certain TF families in Porteresia vis-a- 
vis rice. This unique gene complement might provide 
means of stress adaptation in Porteresia. We identified 
a total of 1 5 1 58 genes responsive to salinity and/or 
submergence stress, which have been implicated in 
diverse biological processes and metabolic pathways. 
The genes involved in similar biological processes, such 
as photosynthesis, secondary metabolite biosynthesis, 
hormone biosynthesis, carbohydrate metabolism, cell 



wall-related enzymes and transport etc., have been 
implicated in salinity and/or submergence stress in 
Arabidopsis, rice and halophytes as vvell.^^'^^'^^''^°''^^ 
A significant number of genes of unknown function 
were also differentially expressed under various stress 
conditions analysed, and many of these genes were pre- 
sumably specific to Porteresia, which might also play a 
crucial role in stress adaptation. 

Several studies have demonstrated that TFs act as key 
regulators of plant response to environmental 
stimuli."^^""^^ The role of certain TF families, such as 
AP2-EREBP, MYB and NAC, in stress responses has been 
well established. We also observed the differential expres- 
sion of members of these TF families under salinity and 
submergence stresses in Porteresia. We detected that 
most of TF genes differentially expressed under salinity 
stress are up-regulated, whereas most of differentially 
expressed TF genes under submergence stress are 
down-regulated. This is in accordance with the total 
number of genes differentially expressed showing up- 
and down-regulation under salinity and submergence 
stresses. These results suggest that overall transcriptional 
activity is induced under salinity stress, but repressed 
under submergence stress. Based on the analyses, we 
found that role of NAC, MYB and WRKY family TFs may 
be more crucial for salinity tolerance, whereas bZIP, 
bHLH, HSF and AP2-EREBP families play a major role 
during submergence stress. The role of SUBMERGENCE1 
(SUB1) QTL in submergence tolerance is well known, 
which harbour a tandem cluster of AP2/ERF family 
genes."^"'"^^'^"'^^ Recently, Niroula etal.^^ proposed the 
existence of SUBl -dependent and -independent 
mechanisms involved in submergence tolerance in wild 
rice plants. However, further investigations are required 
to find the exact mechanism of submergence tolerance 
in Porteresia. 

Metabolic pathways maintain cellular homeostasis 
and are very important for imparting stress tolerance/ 
adaptation in plants. We found significant enrichment 
of genes encoding key enzymes involved in important 
metabolic pathways under stress conditions in 
Porteresia. The production of secondary metabolites 
under stress condition has been well known as adaptive 
feature of plants. We also detected the enrichment of 
biosynthetic genes of certain secondary metabolites, 
such as hydroxycinnamic acid serotonin amides and 
phenylpropanoids etc., under stress conditions in 
Porteresia. Hydroxycinnamic acid serotonin amides 
have only been detected in rice, pepper and a few 
other plants so far and are supposed to act as antioxi- 
dants.^^ Although their exact function is unclear, their 
antioxidant property may help encounter/adapt stress 
condition in Porteresia. Phenylpropanoid biosynthesis 
pathway of initial reactions (cinnamates biosynthesis) 
produces coumarate, which gives rise to a variety of 
secondary compounds such as flavonoids, stilbenes 
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and lignin, which can help in protection against stress 
conditions. 

The plant hormones, ethyleneandABAare well known 
to be involved in adaptive response to environmental 
stresses. Jasmonic acid has primarily been implicated in 
defence response, but is involved in other stress 
responses as well.^'^ Furthermore, the biosynthesis of 
several secondary metabolites by jasmonates has been 
well documented,^^"^^ suggesting that biosynthesis of 
jasmonates may impart stress tolerance in plants via 
transcriptional activation of biosynthetic pathway of im- 
portant secondary metabolites. The enzymes involved 
in cyclic methionine salvage pathway, which ensures 
ethylene biosynthesis via methionine recycling, were 
enriched under salinity stress.^^ We observed up-regula- 
tion of ethylene biosynthesis genes under salinity and 
submergence stresses. Earlier reports have also shown 
the altered ethylene levels underwater, which activates 
genes associated with stress acclimation.'^' ABA is a 
well-known stress implicated in plantstress responses, in- 
cluding salinity stress.^°~^^ The increased accumulation 
of ABA under salinity stress has been reported, which, in 
turn, activates several other genes to provide tolerance. 
However, under submergence stress, the levels of ABA 
are reduced. We also observed the down-regulation 
of ABA biosynthesis genes under submergence stress. 
Recently, it has been proposed that SUB1A gene delays 
leaf senescence via hormonal regulation, which can con- 
tribute to enhanced stress tolerance in rice.^^ 

The enrichment of various pathways related to carbo- 
hydrate metabolism, such as starch degradation, anaer- 
obic respiration and trehalose biosynthesis, is in 
concordance with earlier reports, demonstrating that 
concentrations of soluble sugars increases under 
various abiotic stress conditions.^^'^^ The differential ex- 
pression of genes related to carbohydrate metabolism in 
response to various abiotic stresses has been reported 
previously also.^^ The sugars, sucrose and glucose 
either act as substrates for cellular respiration or act as 
osmolytes to maintain cell homeostasis.^^'^^ In fact, 
one of the present hypotheses for survival of rice under 
submergence stress is based on carbohydrate metabol- 
ism only."^''^''^" Trehalose has been proposed as 
osmoprotectant against various stresses. It has been 
demonstrated that accumulation of trehalose in rice 
plants results in increased soluble sugar content, which 
confers tolerance to different abiotic stresses.^' '^^ 
Furthermore, the elevated transcript levels of the genes 
associated with anaerobic respiration can allow plant 
cellstoinvoketheenergyconserving mode of ATP produc- 
tion undersubmergence stress in Porteres/c/ as speculated 
earlier.' ^'^^ Recently, Miro and Ismail'^ reviewed anaer- 
obic respiration and starch degradation as the main 
traits associated with flooding tolerance in rice genotypes. 

The biosynthetic pathways for suberin and cellulose, 
majorconstituents of cell wall, were also found enriched 



in stress-responsive genes in Porteresia. Although the 
exact mechanism is not known, several studies have sug- 
gested a crucial role of cell wall in stress tolerance.^^ The 
mutations in cellulose synthase genes have also been 
shown to alter stress tolerance in plants.^'^'''^ Suberin is 
a primary constituent of primary cell wall of under- 
ground tissues and is analogous to cutin, which forms a 
layer on plant aerial parts to prevent the direct contact 
of epidermal cell with environment and thus, acts as a 
physical barrier to enhance plant resistance to abiotic 
stress. Based on the analysis of cutin biosynthesis 
mutants, it has been demonstrated that plant cuticle is 
required for the regulation of ABA biosynthesis and 
osmotic stress tolerance.''^ Considering the similar bio- 
chemical nature of suberin and cutin, the increased bio- 
synthesis of suberin may be crucial for stress tolerance in 
Porteresia. Furthermore,thecharacterizationofen/7flncec/ 
submergence 1 (esb1) mutant provided evidence for the 
role of suberin in water balance, ion translocation and 
stress tolerance.^^'^^ Thus, alteration in the degree of su- 
berization may result in improved stress tolerance in 
plants. Overall, it can be concluded that a we 1 1 -organized 
transcriptional regulation is responsible for salinity and 
submergence tolerance in Porteresia. It will be very 
tempting to investigate the precise roleof genes involved 
in the metabolic pathways enriched during salinity and 
submergence stresses in Porteresia. The genetic engin- 
eering of such candidate genes/pathways may confer 
stress tolerance in crop plants. 

In conclusion, we report the transcriptome sequence 
of wild halophyte rice, Porteresia. The results suggest 
that extensive transcriptional reprogramming under 
salinity and submergence stresses is responsible for 
tolerance to these stresses in Porteresia. The transcrip- 
tome analyses presented in this study provide several 
candidate genes involved in salinity and submergence 
tolerance, which can be used as potential targets 
for genetic engineering and opens new avenues for 
future investigations. The functional analyses of genes/ 
pathways identified here along with investigations on 
non-coding regulatory RNAs, epigenetic modifications, 
biochemical complexity and metabolite profiling would 
be very useful to help understand the exact mechanism 
of salinity and submergence tolerance in Porteresia. 

Acknowledgements: We are thankful to the 
anonymous reviewers for their invaluable suggestions. 

Supplementary data: Supplementary data are 
available at www.dnaresearch.oxfordjournals.org. 

Funding 

This work was financiallysupported by the core grant 
from NIPGR. R.G. acknowledges INSPIRE Faculty Award 
from the Department of Science and Technology and 



82 



Transcriptome Analyses in Wild Halophyte Porteresia 



[Vol. 2 1 , 



Innovative Young Biotech nologists Award from the 
Department of Biotechnology, Government of India. 
R.S. acknowledges the Department of Biotechnology, 
Government of India, for research fellowship. 



References 

1 . Walia, H., Wilson, C, Zeng, L, Ismail, A.M., Condamine, P. 
and Close, T.J. 2007, Genome-wide transcriptional ana- 
lysis of salinity stressed japonica and indica rice genotypes 
during panicle initiation stage, Plant Mol. Biol., 63, 
609-23. 

2. Kumari, S., Sabharwal, VP., Kushwaha, H.R., Sopory, S.K., 
Singla-Pareek, S.L. and Pareek, A. 2009, Transcriptome 
mapforseedlingstagespecific salinity stress response indi- 
cates a specific set of genes as candidate for saline tolerance 
in Oryza sativa L, Funct. Integr. Genomics, 9, 1 09-23. 

3. Cotsaftis.O., Plett,D., Johnson, A.A.,etal. 201 1, Root-spe- 
cific transcript profiling of contrasting rice genotypes in 
response to salinity stress, A/lo/. Plant, 4, 2 5-41 . 

4. Mizuno, H., Kawahara, Y, Sakai, H., et al. 201 0, Massive 
parallel sequencing of mRNA in identification of unanno- 
tated salinity stress-inducible transcripts in rice (Oryza 
sativa L.), BMC Genomics, 11,683. 

5. Colmer,T.D. and Flowers, T.J. 2008, Flooding tolerance in 
halophytes. New PhytoL, 1 79, 964-74. 

6. Flowers, T.J. and Colmer, TD. 2008, Salinity tolerance in 
halophytes, New PhytoL, 1 79, 945-63. 

7. Amtmann, A., Bohnert, H.J. and Bressan, R.A. 2 005, Abiotic 
stress and plant genome evolution. Search for new 
mode\s, Plant Physiol., 138, 1 27-30. 

8. Sengupta, S. and Majumder, A.L. 2010, Porteresia coarc- 
tata (Roxb.) Tateoka: a wild rice a potential model for 
studying salt-stress biology in rice, Plant Cell Environ., 33, 
526-42. 

9. Bal,A.R.and Dutt,S.i<. 1 986,Mechanismsofsalttolerance 
in wild rice {Oryza coarctata Roxb.), Plant Soil, 92, 
399-404. 

1 0. Moons, A., Bauw, C, Prinsen, E., Van Montagu, M. and Van 
der Straeten, D. 1995, Molecular and physiological 
responses to abscisic acid and salts in roots of salt-sensi- 
tive and salt-tolerant Indica rice varieties, Plant Physiol., 
107, 1 77-86. 

1 1 . Flatten, J.D., Egdane, J.A and Ismail, AM. 201 3, Salinity 
tolerance, Na+ exclusion and allele mining of HKT1 ;5 in 
Oryza sativa and O. glaberrima: many sources, many 
genes, one mechanism? BMC Plant Biol., 1 3, 32. 

1 2. Miro, B. and Ismail, A.M. 201 3, Tolerance of anaerobic 
conditions caused by flooding during germination and 
early growth in rice {Oryza sativa L), Front. Plant Sci., 4, 
269. 

1 3. Latha, R., Rubia, L, Bennett, J. and Swaminathan, M.S. 
2004, Allele mining for stress tolerance genes in Oryza 
species and related germplasm, Mol. Biotechnol., 27, 
1 01 -8. 

1 4. lnan,G.,Zhang,Q., Li, R,etal. 2004, Salt cress.Ahalophyte 
and cryophyte Arabidopsis relative model system and its 
applicability to molecular genetic analyses of growth 



and development of extremophiles, Plant Physiol., 135, 
1 71 8-37. 

1 5. Oh, D.H., Dassanayake, M., Haas,J.S.,etal. 201 0, Genome 
structures and halophyte-specific gene expression of the 
extremophile Thellungiella parvula in comparison with 
Thellungiella salsuginea {Thellungiella halophila) and 
Arabidopsis, P/flntP/7i/s/o/., 154, 1 040-52. 

1 6. Wu,H.J.,Zhang,Z.,Wang,J.Y,etal. 201 2, Insights into salt 
tolerance from the genome of Thellungiella salsuginea, 
Proc. Natl. Acad. Sci. USA, 1 09, 1 22 1 9-24. 

1 7. Majee, M., Maitra, S., Dastidar, K.G., et al. 2004, A novel 
salt-tolerant L-myo-inositol-1 -phosphate synthase from 
Porteresia coarctata (Roxb.) Tateoka, a halophytic wild 
rice molecular cloning, bacterial overexpression, charac- 
terization, and functional introgression into tobacco- 
conferring salt tolerance phenotype,j. Biol. Chem., 279, 
28539-52. 

1 8. Sengupta, S., Patra, B., Ray, S. and Majumder, A.L. 2008, 

Inositol methyl transferase from a halophytic wild rice, 
Porteresia coarctata Roxb. (Tateoka) regulation of pinitol 
synthesis under abiotic stress. Plant Cell Environ., 31, 
1442-59. 

19. Morozova, O., Hirst, M. and Marra, M.A. 2009, 
Applications of new sequencing technologies for tran- 
scriptome analysis, Anno. Rev. Genomics Human Genet., 
10, 1 35-51. 

20. Jain, M. 201 1 , A next generation approach to the charac- 
terization of a non-model plant transcriptome, Curr. Sc/., 
101, 1435-9. 

2 1 . Jhanwar, S., Priya, P, Garg, R., Parida, S.K., Tyagi, A.K. and 
Jain, M. 201 2, Transcriptome sequencing of wild chick- 
pea as a rich resource for marker development. Plant 
Biotechnol. 10, 690-702. 

22. Patel, R.K. and Jain, M. 201 2, NGS QC Toolkit: a toolkit for 
quality control of next generation sequencing data, PLoS 
One, 7,e3061 9. 

2 3. Garg, R., Patel, R.K., Jhanwar, S., et al. 2011, Gene discov- 

ery and tissue-specific transcriptome ana lysis in chickpea 
with massively parallel pyrosequencing and web resource 
development, P/flntP/!i/s/o/., 156, 1 661 -78. 

24. Kawahara, Y, de la Bastide, M., Hamilton, J. P, et al. 201 3, 
Improvement of the Oryza sativa Nipponbare reference 
genome using next generation sequence and optical 
map data, Rice, 6, 4. 

25. Anders, S. and Huber, W. 2010, Differential expression 
ana lysis for sequence count data. Genome Biol., 1 1, R1 06. 

26. Maere, S., Heymans, K. and Kuiper, M. 2005, BiNGO: 
a Cytoscape plugin to assess overrepresentation of 
gene ontology categories in biological networks, 
Bioinformatics,2\, 3448-9. 

27. Jaiswal, R, Ni, J., Yap, I., et al. 2006, Gramene a bird's eye 
view of cereal genomes, NucleicAcids Res., 34,0717-723. 

28. Garg, R., Sahoo, A, Tyagi, A.K. and Jain, M. 2010, 
Validation of internal control genes for quantitative 
gene expression studies in chickpea {Cicer arietinum L.), 
Biochem. Biophys. Res. Commun., 396, 283-8. 

29. Jain, M., Nijhawan, A., Tyagi, A.K. and Khurana,J.P 2006, 
Validation of housekeeping genes as internal control for 
studying gene expression in rice by quantitative real- 
time PGR, Biochem. Biophys. Res. Commun., 345, 646-5 1 . 



No. 1] 



R. Garg etal. 



83 



30. Agarwal,G.,Jhanwar,S., Priya, P., etal. 201 2, Comparative 
analysis of kabuli chickpea transcriptome with desi and 
wild chickpea provides a rich resource for development 
of functional markers, PLoSOne, 7, e52443. 

31. Ozsolak, F. and Milos, RM. 2011, RNA sequencing 
advances, challenges and opportunities, Nat. Rev. Cenet, 
12,87-98. 

32. Jain, M. 201 2, Next-generation sequencing technologies 
for gene expression profiling in plants, Brief Fund. 
Genomics, 1 1, 63-70. 

33. Gong, Q., Li, R, Ma, S., Indu Rupassara,S. and Bohnert, H.J. 
2005, Salinity stress adaptation competence in the ex- 
tremophile TheUungieiia haiopiiila in comparison with 
its re\at]veArabidopsistlialiana,PlantJ.,44, 826-39. 

34. Munns, R. and Tester, M. 2008, Mechanisms of salinity 
tolerance,/4«nw. Rev. Plant Biol., 59, 651 -81 . 

35. Lee, S.C., Mustroph, A., Sasidharan, R., et al. 2011, 
Molecularcharacterization of the submergence response 
of the Arabidopsis thaliana ecotype Columbia, New 
P/7y to/., 190,457-71. 

36. Voesenek, L.A., Colmer, T.D., Pierik, R., Millenaar, RR and 
Peeters, A.J. 2006, How plants cope with complete sub- 
mergence. New Phytol., 1 70, 213-26. 

37. Vashisht, D., Hesselink, A., Rierik, R., etal. 201 1, Natural 
variation of submergence tolerance among Ara/;/dops/s 
thaliana accessions. New Phytol., 190, 299-31 0. 

38. Wong, C.E., Li, Y, Labbe, A., et al. 2006, Transcriptional 
profiling implicates novel interactions between abiotic 
stress and hormonal responses in Thellungiella, a close 
relative of Arabidopsis, Plant Physiol., 140, 1437-50. 

39. Mustroph, A., Lee, S.C., Oosumi, T, et al. 201 0, Cross- 
kingdom comparison of transcriptomic adjustments to 
low-oxygen stress highlights conserved and plant-specif- 
ic responses. Plant Physiol., 1 52, 1 484-500. 

40. Jung, K.H., Seo, Y.S., Walia, H., et al. 201 0, The submer- 
gence tolerance regulatorSubl Amediates stress-respon- 
sive expression of AP2/ERF transcription factors. Plant 
Physiol.,\^2, 1 674-92. 

41 . Fukao,T,Xu, K., Ronald, PC. and Bailey-Serres, J. 2006, A 
variable cluster of ethylene response factor-like genes 
regulates metabolic and developmental acclimation 
responses to submergence in rice. Plant Cell, 18, 
2021 -34. 

42. Lasanthi-Kudahettige, R., Magneschi, L., Loreti, E., et al. 
2007, Transcript profiling of the anoxic rice coleoptile. 
Plant Physiol. , 1 44, 2 1 8 - 3 1 . 

43. Huang, J., Takano, T. and Akita, S. 2000, Expression of 
alpha-expansin genes in young seedlings of rice {Oryza 
sativa L.),P/of7to, 211,467-73. 

44. Sun,W.,Xu,X.,Zhu, H.,etal. 201 0, Comparative transcrip- 
tomic profiling of a salt-tolerant wild tomato species and 
a salt-sensitive tomato cultivar. Plant Cell Physiol., 51, 
997-1 006. 

45. Sakai, H., Ikawa, H.,Tanaka,T.,etal. 201 1 , Distinct evolu- 
tionary patterns of Oryza glaberrima deciphered by 
genome sequencing and comparative analysis. Plant J., 
66, 796-805. 

46. Chen, J., Huang, Q., Gao, D., et al. 201 3, Whole-genome 
sequencing of Oryza brachyantha reveals mechanisms 
underlying Oryza genome evolution, Nat Commun., 4, 
1 595. 



47. Yamaguchi-Shinozaki, K. and Shinozaki, K. 2006, 
Transcriptional regulatory networks in cellular responses 
and tolerance to dehydration and cold stresses,Annu. Rev. 
Plant Biol., 57, 781 -803. 

48. Nakashima, K., Ito, Y. and Yamaguchi-Shinozaki, K. 2009, 
Transcriptional regulatory networks in response to 
abiotic stresses in Arabidopsis and grasses. Plant Physiol., 
149,88-95. 

49. Bhattacharjee, A. and Jain, M. 201 3, Transcription factor 
mediated abiotic stress signaling in rice. Plant Stress, 7, 
1 6-25. 

50. Fukao,T., Harris, T. and Bailey-Serres,J. 2009, Evolutionary 
analysis of theSu/; 7 gene cluster that conferssubmergence 
tolerance to domesticated rice, Ann. Bot, 103, 143-50. 

51 . Xu, l<.,Xu,X., Fukao,T,et al. 2006, Subl A is an ethylene- 
response-factor-like gene that confers submergence tol- 
erance to rice. Nature, 442, 705-8. 

52. Niroula, R.K., Pucciariello, C, Ho, V.T, Novi, G., Fukao, T. 
and Re rata, P 201 2, SUB1 A-dependent and -independ- 
ent mechanisms are involved in the flooding tolerance 
of wild rice species. Plant]., 72, 282-93. 

53. Jang, S.M., Ishihara, A. and Back, K. 2004, Production of 
coumaroylserotonin and feruloylserotonin in transgenic 
rice expressing pepper hydroxycinnamoyl-coenzyme 
aserotonin N-(hydroxycinnamoyl) transferase. Plant 
Physiol., MB, 346-56. 

54. Wasternack, C. 2007, Jasmonates an update on biosyn- 
thesis, signal transduction and action in plant stress re- 
sponse, growth and development, Ann. Bot, 100, 
681 -97. 

55. Goossens,A., Hakkinen,S.T., Laakso, I., etal. 2003, Afunc- 
tional genomics approach toward the understanding of 
secondary metabolism in plant cells, Proc. Natl Acad. 
Sci. USA, 100, 8595-600. 

56. van der Fits, L. and Memelink, J. 2000, ORCA3, a 
jasmonate-responsive transcriptional regulator of plant 
primary and secondary metabolism. Science, 289, 
295-7. 

57. Sasaki-Sekimoto, Y, Taki, N., Obayashi, T, et al. 2005, 
Coordinated activation of metabolic pathways forantiox- 
idants and defence compounds by jasmonates and their 
roles in stress tolerance in Arabidopsis, Plant J., 44, 
653-68. 

58. Burstenbinder, K., Rzewuski, G., Wirtz, M., Hell, R. and 
Sauter, M. 2007, The role of methionine recycling for 
ethylene synthesis in Arabidopsis, Plant]., 49, 238-49. 

59. Fukao, T, Ye ung, E. and Bailey-Serres, J. 201 2, The sub- 
mergence tolerance gene SUB1A delays leaf senescence 
under prolonged darkness through hormonal regulation 
in nee. Plant Physiol., 160, 1 795-807. 

60. Xiong, L. and Zhu, J.K. 2003, Regulation of abscisic acid 
biosynthesis, P/flnt P/ti/s/o/., 133, 29-36. 

61 . Verslues, RE. and Zhu, J.K. 2005, Before and beyond ABA 
upstream sensing and internal signals that determine 
ABA accumulation and response under abiotic stress, 
Biochem. Soc. Trans., 33, 3 75-9. 

62. Tuteja, N. 2 00 7, Mechanisms of high salinity tolerance in 
plants. Methods £«zi/mo/., 428, 41 9-38. 

63. Cutler, S.R., Rodriguez, PL, Finkelstein, R.R. and 
Abrams, S.R. 201 0, Abscisic acid emergenceof a core sig- 
naling network,An«u. Rev. Plant Biol., 61, 651 -79. 



84 



Transcriptome Analyses in Wild Halophyte Porteresia 



[Vol. 2 1 , 



64. Sail<a,H.,Okamoto,M.,Miyoshi,l<.,etal. 2 00 7, Ethylene pro- 
motes submergence-induced expression of OsABAS ox 1 , a 
gene that encodes ABA 8'-hydroxylase in rice, Plant Cell 
Physiol., 4, 287-98. 

65. DubeyR.S. and Singh, A.K. 1 999, Salinity inducesaccumu- 
lation of soluble sugars and alters the activity of sugar me- 
tabolizing enzymes in rice plants, Biol. Plant, 42, 233-9. 

66. Gill, P.K., Sharma, A.D., Singh, P. and Bhullar, S.S. 2001 , 
Effect of various abiotic stresses on the growth, soluble 
sugars and water relations of sorghum seedlings grown 
in light and darkness, Bulg.j. Plant Physiol., 27 , 72-8^. 

67. Seki, M., Narusaka, M., Ishida, J., et al. 2002, Monitoring 
the expression profiles of 7000 Arabidopsis genes under 
drought, cold and high-salinity stresses using a full- 
length cDNA microarray Plant]., 31,279-92. 

68. Gupta, A.K. and Kaur, N. 2005,Sugarsignallingand gene 
expression in relation tocarbohydrate metabolism under 
abiotic stresses in plants,/. Biosci., 30, 761 -76. 

69. Rosa, M., Prado, C, Podazza, G., et al. 2009, Soluble 
sugars— metabolism, sensing and abiotic stress a complex 
network i n the I ife of plants, PlantSignal. Behav. , 4, 3 8 8 -9 3 . 

70. Voesenek, L.A. and Bailey-Serres, J. 2009, Plant biology 
genetics of high-rise rice. Nature, 460, 959-60. 

71 . Garg, A.K., Kim, J. K., Owens, T.G.,etal. 2002, Trehalose ac- 
cumulation in rice plants confers high tolerance levels to 
different abiotic stresses, Proc. Natl. Acad. Sci. USA, 99, 
1 5898-903. 



72. Redillas, M.C.F.R., Park, S-H., Lee, J.W., et al. 2012, 
Accumulation of trehalose increases soluble sugar con- 
tents in rice plants conferring tolerance to drought and 
salt stress, Plant Biotechnol. Rep., 6, 89-96. 

73. Koiwa, H. 2009, Pathways and genetic determinants for 
cell wall-based osmotic stress tolerance in the Arabidopsis 
thaliana root system. In: Jenks, M.A. and Wood, A.J., 
editors. Cenes for plant abiotic stress, Ch. 2. Wiley- 
Blackwell, Oxford, UK. 

74. Chen,Z., Hong,X.,Zhang, H.,et al. 2005, Disruption of the 
cellulose synthase gene, AtCesA8/IRX1 , enhances 
drought and osmotic stress tolerance in Arabidopsis, 
P/flnt;., 43, 273-83. 

75. Zhu,J.,Lee,B.H.,Dellinger,M.,etal. 201 0,Acellulosesyn- 
thase-like protein is required forosmotic stress tolerance 
in Arabidopsis, Plant]., 63, 1 28-40. 

76. Wang, Z.Y., Xiong, L, Li, W, Zhu, J.K. and Zhu, J. 20 1 1 , The 
plant cuticle is required for osmotic stress regulation of 
abscisic acid biosynthesis and osmotic stress tolerance 
in Arabidopsis, P/fl«/:Ce//, 23, 1971 -84. 

77. Baxter, 1., Hosmani, PS., Rus, A., et al. 2009, Root suberin 
forms an extracellular barrier that affects water relations 
and mineral nutrition in Arabidopsis, PLoS Genet., 5, 
el 000492. 

78. Franke, R.B., Dombrink, I. and Schreiber, L. 201 2, Suberin 
goes genomics use of a short living plant to investigate a 
long lasting polymer, Front. PlantSci., 3, 4. 



